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We present a general scheme for treating particle beams as many particle systems. This includes 
the full counting statistics and the requirements of Bose/Fermi symmetry. In the stationary limit, 
i.e., for longer and longer beams, the total particle number diverges, and a description in Fock space 
is no longer possible. We therefore extend the formalism to include stationary beams. These beams 
exhibit a well-defined "local" counting statistics, by which we mean the full counting statistics of all 
clicks falling into any given finite interval. We treat in detail a model of a source, creating particles 
in a fixed state, which then evolve under the free time evolution, and we determine the resulting 
stationary beam in the far field. In comparison to the one-particle picture we obtain a correction 
due to Bose/Fermi statistics, which depends on the emission rate. We also consider plane waves 
as stationary many particle states, and determine the distribution of intervals between successive 
clicks in such a beam. 
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I. INTRODUCTION 

In the standard framework of quantum mechanics each 
type of systems is assigned a Hilbert space H. Prepara- 
tions are given by density operators on H, and measure- 
ments are given by operator valued measures on some 
outcome set. The basic predictions of the theory are 
probabilities, i.e., the asymptotic relative frequency, to 
be determined in a sufficiently long run experiments with 
the same preparation and measurement. However, many 
experiments do not follow this standard scheme: instead 
one often prepares a beam of particles. On the mea- 
surement side one then uses detectors, whose readout is 
a time series of clicks. Instead of probabilities the ba- 
sic experimental results are then count rates. There are 
many reasons to study the relationship between these two 
models of quantum experiments. We were led to it in the 
course of a project on the tomography of single photon 
sources. Clearly, in this case it is relatively easy to get 
beam-type data, in which, for example, the antibunching 
dip in the -correlation function indicates single parti- 
cle emission [IJ[2]. On the other hand, the emission time 
of individual particles is not controlled in beam mode, 
so absolute click times are meaningless. Hence it is im- 
possible to get a full tomography of the single particle 
state created by such a source on demand; beam data 
and single-shot data have to be combined. 

It is clear that beams can be described as many-particle 
quantum systems, and it is natural to use the standard 
Fock space for Bosons and Fermions. This will indeed be 
our starting point. However, the Fock space arena is too 
restrictive for describing truly stationary beams, which 
necessarily contain infinitely many particles. The full 
counting statistics of stationary beams will be obtained 
by taking the limit of longer and longer finite beams in 
a suitable way. The resulting description will be easier 
than that of finite beams, which must implicitly always 
contain the particulars of switching the beam on and off. 
Demonstrating this, and providing a theoretical frame- 
work for stationary beams is the main goal of this pa- 



per. On the measurement side this involves observables, 
whose outcome space is a (usually infinite) set of count- 
ing events. They are constructed directly in terms of one- 
particle observables, and do not rely on the identification 
of field intensities and count rates made in quantum op- 
tics via the Glauber model [3] . On the side of the sources 
we focus on the type of sources, which make the connec- 
tion to the one-particle picture most apparent, and which 
best realize the idea of a train of independent particles. 
Technically, these are gauge invariant quasi-free states. 
We are, of course, aware that that this excludes many 
interesting phenomena, particularly for photons. But in 
our ongoing research on the subject we found the quasi- 
free beams to be an important starting point for more 
complex situations, e.g. by conditioning. 

In order to make the paper more self-contained, we in- 
clude brief statements of the relevant prerequisites, like 
the theories of counting processes 0HB], Fredholm de- 
terminants 7 , and quasi- free states [H [§]■ Processes 
rather similar to the ones we find have been discussed 
in the mathematical literature under the heading of "de- 
terminantal processes" [TU] . However, to the best of our 
knowledge, these were always considered to be station- 
ary in space I12j rather than in time. Stationarity 
in time requires some extra work, e.g., bringing in co- 
variant arrival time observables on the one particle level, 
and the second quantization of generalized (POVM) ob- 
servables. In the end, an appropriate "local trace class" 
condition characterizing the particle sources can be for- 
mulated quite simply also in the time case. In this we 
see the main contribution of our paper, although we also 
hope that for some readers it will serve as an invitation to 
a more systematic view of quantum counting processes. 

Our paper is organized as follows. In Sect. |Hj we in- 
troduce some notation, explain the concept of a point 
process, describe the construction of general counting ob- 
servables, and give an account of arrival time measure- 



ments. In Sect. Ill we introduce the class of quasi-free 
states in Fock space, and compute the full counting statis- 
tics for such states. We also give the essential technical 
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result which makes it possible to extend the results to 
the local statistics for stationary beams. Sect. IV is de- 



voted to a concrete way of going to the stationary limit 
starting from an explicit dynamical description of parti- 
cle creation. In Sect. [V] we give a general scheme for de- 
scribing local counting statistics for quasi-free stationary 
beams; in particular, we get explicit formulas for correla- 
tion functions. Finally, in Sect. |VI[ we apply the results 
to the simplest example of a stationary beam, the plane 
wave viewed as a many-particle state. Often enough such 
an interpretation is suggested in textbook treatments of 
scattering solutions of the stationary Schrodinger equa- 
tion for one-particle potential scattering. Here we take 
it literally, and as a bonus get correlations and waiting 
time distributions in such a beam. 



II. COUNTING OBSERVABLES 

Let us fix some notation. Throughout, % will be the 
Hilbert space of a single particle. By T S (H) we denote 
the Fock space over H, with s = +1 for Bosons and 
s = — 1 for Fermions. When it is irrelevant, or clear from 
the context, the index s will be omitted. For an operator 
A on T~L, we will denote by r s (A) the operator, which 
on TV-particle wave functions acts like the TV-fold tensor 
power A® N . Clearly, T S {AB) = T S (A)T S (B). 

For operator valued or scalar valued measures \x we 
abbreviate the integral over a scalar function / as 



(i) 



i.e., we use round brackets for the set function and square 
brackets for the integral, i.e., the expectation value func- 
tional in the case of a probability measure. 

For observables (POVMs= "positive operator valued 
measures") the discrete case, in which points have fi- 
nite measure, is often used. Then F[f] = ^ x f{x)F x , 
where F x are positive operators with F[\] — F% — 1- 
The projection valued special case is characterized by 
F x F y — 5 xy F x , or in a form also valid in the continuous 
case: F[fg] = F[f]F[g]. 

In most textbooks, observables are simply identified 
with self-adjoint operators A, which presupposes that 
X C K and takes F as the spectral measure of A, so 
A = jF(dx)x — F[x]. The same measure also defines 
the one-parameter unitary group exp(itA) = F[e ltx ~\ 
generated by A. Generators are second quantized by 
A i ^ dT(A), where T(exp(itA)) = exp(itdT(A)) for all 
t. However, for the purpose of this paper it is much more 
appropriate to start from POVMs on some outcome space 
X, which need not be the real line. The point is that for 
the natural second quantization of such POVMs [T3] the 
outcome space also changes: where the observable F at 
the one-particle level gives the probability of outcomes 
x £ V, its second quantization TF will correspond to a 
measurement of F on every particle, and the result of 
this measurement is a distribution of points in X. The 



probabilities for such measurement outcomes constitute 
a so-called "point process". 



A. Point processes 

A point process is a probability measure on the space 
of outcomes, where each outcome is a collection of (not 
necessarily distinct) points in a set X. We can think 
of each outcome as a possibly infinite numbered list 
(x(l),x(2), ...) with x{i) £ X, with the understanding 
that the ordering of the elements is irrelevant but, in 
contrast to the set {x(l), x(2), ...}, we do count the num- 
ber of occurrences of each x £ X. A way to express this 
compactly is to take as the outcome of a point process 
the so-called empirical measure 



x(i) ) 



where S x denotes the point measure (with (5-function den- 
sity) at x £ X. A measure £ of this form is also called 
a counting measure, and is characterized among mea- 
sures by the property that the measure of each set is an 
integer, namely the number of points Xi in that set. Us- 
ing the bracket notation introduced above, we then have 
£[/] = J2i f( x ("i)) for an empirical measure £ = £\ S x{i) . 
Since this bracket is linear in /, we can use it to charac- 
terize the probability distribution of a point process by 
its Fourier transform, i.e., by the expectation (•) of the 
function £ i— >• exp(i£[/]). This is called the characteristic 
function 



C(J) 



(2) 



of the distribution and contains the full counting statis- 
tics. For example, consider k disjoint subsets Xp of X, 
and let \p denote the indicator function of Xp. Then 
we get the characteristic function of the joint probability 
distribution for the number of counts in the sets Xp as 



E 



p(m, . ..n k )e 



= C(£,A, X£ ). (3) 



If the particle numbers are independent for every parti- 
tion of X into sets Xp, we have a Poisson process, which 
is characterized by a measure (i on X, called the in- 
tensity measure of the process, such that p(ni, . . . n^) — 
Y\ e ji(Xp) n /(nl) exp(—fi(Xp) and hence, 



C(/)=exp L{dx){e^ -I). 



(4) 



The k th moment of the point process is defined as the 
uniquely determined permutation symmetric measure mk 
on X k , satisfying 



rrik(dxi 



■ dx k ) JJ f(a 

3=1 
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Using arbitrary functions / on X k , we can equivalently 
give the definition as 



m k[f] = ( E /(^l)' ' ' • > X fe)) 



(5) 



Since 



c{f) = E 



fc! 



m k {dxi ■ ■ ■ dx k ) Y\_ f( x j)> ( 6 ) 

3=1 



we can extract the moments from C(A/) by differenti- 
ating with respect to A. For a Poisson process, the ex- 
pansion of the characteristic function in powers of f is 
C(f) = */*[/] - \{n[P] + M[/f ) + 0(/ 3 ), so the first mo- 
ment is mi = \x and the second is m-i = [i ® fi + /i o A -1 
with the diagonal map A(x) = (x,x). The second mo- 
ment thus has a singular part concentrated on the diag- 
onal. This is not a special feature of the Poisson process, 
but occurs for any counting process. It is therefore cus- 
tomary to consider a modified set of moments, called 
factorial moments [2 E] , or "correlation functions" [10] , 
which do not have such singularities. Like the moment 
TOfc, the factorial moment of order k, which we denote by 
fhk, is a permutation symmetric measure on X k . For a 
function / of k variables, the factorial moment is defined 
by the following expectation: 



m k 



;[/] = ( E f{x{ii),...,x{ik)) 



(7) 



distinct 



By comparing the expression Q to ([5j it is clear that 
the exclusion of multiply occurring indices in the former 
just has the effect of eliminating the singular term from 
the second moment. For the Poisson process one has 
fh k = n® k for all k. 

The factorial moments are most easily obtained from 
the characteristic function C by observing that the gen- 
erating function 

k 

e w = e h I ^(^i • • • dx ^ n f( x i) ( g ) 

k ' J 3=1 

is related to C just by a transformation of the argument: 
C(f) = Cie** - 1) (9) 

As an example of using ([9]), consider the probability 
PY{n) of finding exactly n particles in a measurable re- 
gion Y G X. By ([3| and we get the relation 



(10) 



for z = e lA , where \y is the indicator function of Y. 
In particular, we get the no event probability directly 



from the factorial moment generating function by ana- 
lytic continuation: 



Py(0) = C(-xy) 



(11) 



These probabilities determine the interval statistics of 
a point process on the time axis X = K.. Indeed, let 
Po(tii — P[ti,t 2 )(0) denote the probability for not find- 
ing a click in the interval [ti,ta). Then the probability 
of having no click on [fi,f2) and at least one click just 
before t\, say in the interval \t\ — s,ti), is the same as 
having no click on [£i,t2) and at least one in [ii — e, tz), 
which is po(*ij ^2) — Po{ti — e, ^2)- Hence the conditional 
probability for having no click on [t\, t%), on the condition 
of having at least one in [t\ — e,ii), is 

PojtiM) -Po{h - e,t 2 ) 
1 -po(h - e,ti) 

At the limit e — > 0, this tends to 1— the probability of 
having to wait at most time r = ti — 1\ for the next click, 
after a click at t±. Hence the probability density w s (t) 
for the waiting time r 6 [0, 00) given a click at s is 



w s (t) 



dp 
dt x 



d 2 Pa 
dhdt 2 



(8,8 + T). (12) 



Since we will eventually apply the above formalism to 
particle detection processes, we close this subsection by 
a remark on the role of the factorial moment densities 
in the standard theory of photon counting used in quan- 
tum optics. If X has a natural measure dx (typically 
X = R m with the Lebesgue measure) , we can often write 
fhk{dx\ . . . dxk) — h^ k \x\, . . . ,Xk)dx\ . . . dxk for a den- 
sity function h^ k \ In Glauber's model of a photon detec- 
tion process, we have 



t ) = const • G^ n '(xi, ...,x n 



,xi), 



where is the usual "correlation function" defined 

using the field operators QJQ3]. I n this context, one also 
typically uses the normalized correlation functions, which 
we define for a general point process by 



g {n) {x ll ...,x n ) 



h^(x u ...,x n ) 

hW(xi)>--hM(x n y 



(13) 



B. Second quantization of general observables 

One reason for introducing characteristic functions of 
point processes is that they make the construction of the 
second quantized observable TF from the single particle 
observable F extremely simple |13j . Indeed, if we just 
express the idea that TF measures F on all the particles 
we get, restricted to the iV-particle space, the operator 



((TF) 



Mf] 



A'* 



F(d Xl ) ® ••• ® F(dx N )e 



= (F[e*f\) 
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Taking the direct sum over N, we get the fundamental 
formula 



(TF) 



(14) 



This expression makes sense on full Fock space, i.e., it 
does not require restriction to the Bose or Fermi sector. 
Hence, for a state given by a density operator p on this 
space the full counting statistics can be extracted from 
the characteristic function 



C(f)=trp(TF) 



trpr(F[. 



(15) 



Of course, for a Bose or Fermi system, the operator p has 
support in the appropriate subspace and we can replace 
the r in this formula by the corresponding restriction T s . 

The factorial moments can be computed from this us- 
ing Eqs. Q and For a state p on full Fock space we 
get 

C(/)=tr P r(l + F[/]) 

oo 

= Y,trp N (l + F\f\)* N 



N=0 

oo 

E 

fc=0 



trp k F[ff 



(16) 



where pk is the reduced fc-particle reduced density oper- 
ator 



E 

N=k 



tr 



[k + l,...,N] PN- 



(17) 



Hence the factorial moments are 
fhk{dxi • ■ ■ dxk) = k\ tr pkF(dxi) 



F{dx k ). (18) 



C. Arrival time observables and their dilations 

Due to an old argument of Pauli, an arrival time ob- 
servable cannot be a spectral measure of a self-adjoint 
"time operator" . However, the generalization of the no- 
tion of observables to POVMs immediately allows time- 
shift covariant observables to be constructed [TS] . In this 
subsection we describe the general construction of ob- 
servables, which measure the arrival time t and arrival 
location a: of a particle [16] , Here location is taken in a 
rather broad sense, and could just be the number of the 
detector which responds. We consider arbitrary observ- 
ables, which are covariant for time translations, i.e., 



cxp(iHt)F [/] exp(-iHt) = F[nf], 



(19) 



where H is the Hamiltonian, and r is the time shift on 
functions of t and x, i.e., (r t /)(x,t') = f(x,t' — t). The 
standard method [16] to build all covariant observables 
(even for a general covariance group with representation 
g i y U g ) involves two steps: one first uses the Naimark 



dilation theorem to turn any generalized (POVM) ob- 
servable into a projection valued one, say F, which lives 
on another Hilbert space T~L and is connected to F by 
an isometry V : U -> U so that F[f] = V*F[f]V. 
There is also a unitary group representation U on T-L 1 
for which F is covariant, and which is intertwined by V, 
i.e., UgV — VU g . In the second step one uses the the- 
ory of Mackey [17] who called projection valued covariant 
observables "systems of imprimitivity" and showed their 
intimate connection to induced representations. This sec- 
ond part is easy for just the time translation group R, 
and leads to standard Schrodinger pairs of position and 
momentum operators, with some multiplicity. Thus in 
the dilation space "energy" is the canonical multiplica- 
tion operator canonically conjugated to "time" and has 
therefore purely absolutely continuous spectrum. The co- 
variant time observable approach is therefore limited to 
Hamiltonians H with absolutely continuous spectrum, in 
which case the Hilbert space is of direct integral form 



H 



dEU 



(20) 



This is shorthand for the space of wave functions, which 
are functions of energy such that ip(E) e He, the multi- 
plicity space at E. This will, of course, be {0} when E 
is not in the spectrum of H (e.g. when E < for the 
standard kinetic energy). Scalar products are computed 



as 



<</#)„ = / dE (cb(EM(E)) 



He 



(21) 



with the scalar product of T-Le- The technical (measura- 
bility) conditions on direct integral Hilbert spaces are to 
ensure that this expression makes sense (see e.g. |18j). 
Of course, the Hamiltonian is the multiplication oper- 
ator (Hil))(E) = Etjj(E) in this representation. More 
generally, a bounded operator A commutes with H if 
(Aip)(E) = Ae"4>(E) for some measurable family of oper- 
ators Ae € B(He)- We write for this 



A= / dEA, 



(22) 



If a space (20 1 allows a projection valued covariant 



time observable, and hence a self-adjoint conjugate time 
operator, this operator generates a unitary group which 
shifts the energy variable. It thus introduces a canoni- 
cal identification between all the spaces He- In partic- 
ular, they must be non-zero also for negative energies, 
which is exactly the above-mentioned argument of Pauli 
that semi-bounded Hamiltonians do not allow a projec- 
tion valued time observable. Nevertheless, this structure 
appears as the dilation of any given time observable. The 
Hilbert space in that case can be written either as the 
tensor product £ 2 (K, dE) <g> JC or, in the spirit of (20 1, 
as the space of /C-valued £ 2 -functions on K. The time 
observable in this case is computed in the usual way by 
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Fourier-transforming to C 2 (M., dt)®/C, and the joint mea- 
surement of t and x is realized in this tensor product. It 
is thus characterized by the following data: 

1. a Hilbert space IC, which will be the energy- 
independent multiplicity space of the dilated ob- 
servable, 

2. a family of isometries Ve : He ^ IC, which together 
define the dilation isometry V : U -> £ 2 (R, dE)®K 
via (Vtp)(E) = V E ip(E), and 

3. an observable G with outcome space X in the 
Hilbert space IC. 

We have to compute the expectation operator F[h] for 
arbitrary functions / of (t, x) but it suffices to do this for 
the product functions f(t,x) — h(t)g(x). For these the 
above data determine the operator 

F[hg]1>(E) = fdE' h(E - E') V* E G[g]V E ^(E'), (23) 



where h denotes the Fourier transform of h, normalized 
as 



h(E) = I dt e lEt h{t\. 



(24) 



This ensures that for g = 1 and /i/lwe find F(h) /* I, 
so the observable F is normalized. It is convenient to 
allow also subnormalized observables, i.e., F(l) < 1. In 
that case the operator 1— F(l) measures the probability 
that the particle never arrives. By construction this op- 
erator will commute with H, and the only modification 
in the above setup is to allow V to be a general operator 
with ||V|| < 1, rather than an isometry. 



III. QUASI-FREE STATES 

A. Physical background 

Let us begin with a Boltzmann statistical model of 
a multi-particle preparation: Suppose we have a one- 
particle preparation with density operator <7q, which we 
run at N random times ij. Hence, if a(t) = e~ lHt <7oe lHt 

is the time translate of ao, we get the state &) i=1 <y(U). 
We ignore for the moment the symmetrization require- 
ments, so we apply the observable F to each of these 
systems separately, obtaining a point Xi as a measur- 
ing result. In order to determine the characteristic func- 
tion of the counting statistics we need the distribution of 
the emission times, which we take to be Poisson with 
intensity measure /i (i.e., with characteristic function 
Ctime(ff) = exp J /j,(dt)(e 19 — 1), see Q). From this we 
get the characteristic function of the counts 



C(/)=exp / fi(dt)tva(t)F(e if - 1). 



This is again Poisson, and depends only on the integral 
(7 = f (i(dt)cr(t). The corresponding state on full Fock 
space is Q) N -^ycr® 7 ^, where the factorial is the usual cor- 
rection factor for indistinguishability familiar in classical 
statistical mechanics. 

Of course, this state is not consistent with Bose or 
Fermi statistics. Its closest analogue is to replace -mcr® N , 
by P s a® N P s , where P s denotes the projection onto the 
(anti-)symmetric subspace oiW® N . That is, we consider 
the quasi-free state with density operator 



trT s (a) 



(26) 



on r ' sl-L. Here we do not take <r to be normalized. Instead 
the normalization factor of a determines the particle 
number distribution. To be precise, (26) is a "gauge in- 
variant" quasi-free state. More general quasi-free states, 
which do not necessarily commute with particle num- 
ber are defined in 19J, and will not be studied in this 
paper. Quasi-free states are plausible models for non- 
interacting particle beams, because they give the same 
results as Boltzmannian independence in the weak beam 
limit. They also describe naturally the production of a 
beam. Consider an oven, modeled as an ideal gas with 
one-particle Hamiltonian H, at temperature T = l/(kj3) 
and chemical potential \x. Then we have the grand canon- 
ical quasi-free state with a — exp(— j3(H — fit)). We then 
get a beam by letting some particles escape through a 
hole, and we can also add (possibly time-dependent) one- 
particle potentials, collimating filters and the like. The 
important point is that as long as we only apply one- 
particle operations, i.e., unitary operators of the form 
r([/), the quasi-free character of the initial state will be 
preserved. 

In the context of quantum optics this type of photon 
beam is usually called thermal light (see e.g. [1]), the 
state appearing as a special case of "chaotic state" |14) . 
The latter is defined in terms of the occupation number 
states |{?i/c}) as 



En 

{n k } k 



\{n k })({n k }\, (27) 



where k indexes the modes, and the a k > are parame- 
ters satisfying a k < oo, each a k coinciding with the 
expectation value of the mode k occupation number. In 
fact, any quasi-free state of the form ( 26 ) for s = 1 can 



be written as (27) by choosing the modes according to 



(25) 



an eigenbasis of the positive trace class operator tr; the 
parameters a k are then the eigenvalues of a = a — a). 



B. Characteristic functions 

A crucial tool in the following is a formula for the de- 
nominator in (26). When A is trace class (i.e., \\A\\i = 
tr vMt/1 < oo) 



and, in the Bose case \\A\\ < 1, then 



G 



T S (A) is also trace class and 

trT s (A) = dct(l - sAY 



(28) 



For the theory of such infinite dimensional determinants 
we refer to [7J. Using this formula, we get a simple ex- 
pression for the characteristic function of a counting mea- 
surement: 



C(f) = 



tiT(o-)r(F[e i/ ]) _ trT(aF[e if ]) 



trT{a) 
det (1- saF[e l f]y 



trr(cr) 



det (I - sa)- s 
= det (1 - s(l - sa)~ s a (F[e if ] 

To summarize: 



I))" 



C(f) = det (1 - sdF[e if - l])~ 



with 



1 - sa 1 



(29) 
(30) 



where we took the liberty to write a fraction because nu- 
merator and denominator commute, and the expression 
can be evaluated in the functional calculus. It is useful 
to note the bounds on the operators a, a in the Bose and 
Fermi case: Clearly both operators must be positive and 
have finite trace. In the Bose case we need in addition 
that a < (1 — e)l for some e > 0, which is equivalent 
to saying that a is bounded. In the Fermi case it is the 
other way around: a can be any bounded operator, which 
implies that u is strictly less than the identity. The for- 
mula (29 1 contains the complete counting statistics for 



the counting observable (compare also [2"U]). 



Factorial moments 



From ( 29 ) and ^ we get the factorial moment gener- 
ating function 

C(f) = det (i - sdF[f}y s = trT s (aF[f]) (31) 



Comparing this with ( 16 ) we see that the fc-particle re- 



duced density operators of the quasi-free state are p k — 
P s a® k P s , so by (18 1, the factorial moments are given by 



m k (dxi ■ ■ -dx k ) = k\tv P s a® k P s F(dxi) <g> • • ■ F(dx k )- 

(32) 

The first moment is simply fh\(dx) — tr aF(dx); for 
the second moment, the expression can be further re- 
duced, so that traces have only to be taken in the 
one particle space. To this end we write the (anti- 
)symmetrization projection P s = (1 + sF)/2, where F is 
the unitary transposition operator, and use tr(F A®B) = 
tr(AB). Then 



rhi(dx dy) — fhi(dx)fhi(dy) 

+s ti(aF(dx)aF(dy)). 



Now the trace on the right hand side is a positive mea- 
sure on X x X and is also positive definite in the sense 
that it gives positive expectation to functions of the form 
f(x)f(y). This shows that for Bosons we always have 
the bunching effect g^ > 1 and the antibunching effect 
q(2) < i f or Fermions. (Recall the definition (13) of the 



correlation function g^). Clearly, there are interesting 
cases of photon antibunching, but these require artfully 
correlated, not quasi-free sources. 

For later use we note the fc th order generaliza- 
tion of (33). The expression for tr(V 7r A 1 <g> -A k ), 
for a permutation operator V T is based on the cy- 
cle decomposition of the permutation ir, say ir — 
(«!,..., i r )(ji, ■ ■ ■ , j s ) ■ ■ ■ , and gives the product of the 



traces tr( J 4 il • • • A ir ) tr(Aj 1 ■ ■ ■ Aj s ) 
to introduce the measures 



It is convenient 



Ht(dxi ■ ■ ■ dxi) = tr ( JJ dF(dx a ) 



(34) 



so that Hi(dx) — fhi(dx), and (33) reads fhi(dxdy) — 
fii(dx)fi2(dx)+sfi2(dxdy). Then, for example, for k — 3, 
we find 

fha(dxdydz) = fxi(dx)fxi(dy)fj,i(dz) 

+s (n<2(dx dy) n\(dz) + cyclic) 
+2$iefi 3 (dxdydz). (35) 

For general k we get similar expansions into products 
of measures, the combinatorics of which requires some 
representation theory of the permutation group, which 
we will not expound here. 



D. The localization Lemma 

The main aim of our paper is to establish local count- 
ing statistics even in situations where the global particle 
count is infinite, as will be the case for any stationary 
beam, or translationally invariant gas. The idea is to 
use the formula ( 29 1 for the characteristic function even 



in situations, where the operator a has infinite trace, but 
the product aF is sufficiently well behaved so the formula 
makes sense as written. Actually, even more general sit- 
uations can be covered, if we replace aF by \faF\fa. 
That is we start from the characteristic function 



C(f) 




sVdF[e lf -l]Vd 



(36) 



(33) 



This is the same as ( |29[) w hen a has finite trace. Indeed, 
we have the identitvpT] det (I + AB) = dct(l + BA), 
whenever both AB and BA are both trace class. In the 
case at hand this is applied to the two Hilbert-Schmidt 
operators A = \fa and B — F[e l f — The form 

given in the following Lemma is even slightly more gen- 
eral. 

Lemma 1 Let F be a measure on some set X, whose 
values are positive operators on a Hilbert space Hi, with 
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F(X) < 1, and consider a subset Xq C X . Let 'H2 be an- 
other Hilbert space and W : H2 — ► "Hi a bounded operator 
such that 



tr W*F{X )W < 00, 



(37) 



and in the Fermi case (s = —1) also \\W\\ < 1. Then the 
formula 

C{f) = det (1 - sW*F[e if - l]wy s , (38) 

for all f vanishing outside Xq defines the characteristic 
function of a point process in Xq . 



Xq C X . Let Woo, W a : H2 — > "Hi be bounded operators, 
such that W* — > strongly, and 



w a w*<w w*. 



(40) 



Then Woo and each W a satisfy the conditions of Lemma 
[7J and for the associated characteristic functions C a , Coo 
holds 

C 0O (/)=lim(7 a (/) 

a 

uniformly for f with support in X . 



Proof Consider the Naimark dilation F — V*FV. Then 
for / with support in Xq we can write 

W*F[f]W = W*V*F[f]VW 

= W*V*F(Xo)F[f}F{Xo)VW 
= W*F[f]W, 

where W = F(Xq)VW , and at the second equality we 
used the projection valuedness of F. Now by assump- 
tion (37 1 the operator W*W has finite trace, i.e., W is 



a Hilbcrt-Schmidt operator. Moreover (relevant only for 
s = -1) ||W|| < 1 implies \\W\\ < 1, because F(X ) 
is a projection, and the dilation operator V satisfies 
V*V = F(X) < I. Hence WW* satisfies all conditions 
required of an operator a to define the counting statistics 
of a bona fide quasi-free state, with respect to a second 
quantized observable. The associated characteristic func- 
tion ([291) is 



C(f) = det (l - sWW*F[e lf - l] 



(39) 



which has the form stated in the lemma by the same 



argument that gave the equality of (29) and (36) at the 
beginning if this section. ■ 

Since we can take W = yfa in the Lemma, we find 
tr y/aF(X ) yfa < 00 as a sufficient condition to apply 
(36). The similar looking condition ||ctF(Xo)||i < 00, 



which is suggested by the characteristic function ( 29 1 , 
is actually stronger. This is implied by the estimate 
tr y/AByfA < ||AB||i, which holds for arbitrary positive 
operators A, B. (For a proof note that the trace norm 
is the sum of the singular values, which dominates the 
sum of the absolute values of the eigenvalues [71 Thm. 
1.15], and that AB and y/ ABy/A have the same nonzero 
eigenvalues.) 

As a byproduct we can now approximate the counting 
statistics of a stationary state by the counting statistics 
of finite-beam ones: 

Lemma 2 Suppose that Wq : H2 — ► "Hi satisfies the con- 
ditions of Lemma [7] for a given measure F and a set 



Proo f D efining W a for W a as in Lemma [lj the assump- 
tion (40) gives W a W* < W Wq, which implies that the 
conditions of Lemma [l] are valid also for each W a and 
Woo, and that ||W*^|| 2 < ||W *V>|| 2 for ip £ H where 
now Wq is Hilbert-Schmidt. Together with the strong 
convergence W* — > W^, this implies convergence in the 
Hilbert-Schmidt norm. Hence, W a W* — > WooW^ in the 



trace norm. Using ( 39 1 , we get the uniform convergence 
of the C a from the estimate | det (I + A) - det(l + B)\ < 
\\A-B\\ 1 e^+^ +1 [21]. ■ 

Now, in order to approximate a stationary beam with 
non-trace class a, by a sequence of finite beams with 
trace class a a , it is sufficient to have the weak operator 
convergence yfa a — > y/a, and majorization a a < cto by 
some <7q with tr y/aoF(XQ)y/aQ < 00. This will give uni- 
form convergence of characteristic functions for the local 
counting statistics. 



E. "Small s" and parastatistics 

We end this section with a remark on parastatistics 
and weak beams. If we write ( 29 ) as 



C s {f) = exp tr — log(I - sd?F [e^ - l] ) , (41) 



the formula gives corrected Boltzmann statistics (251 for 



s = 0, and parafermi (resp. parabose) statistics of order 
p for s = —l/p (resp. s = l/p). From (41) we get a 
rather uniform notion of weak beams: Whenever a is 
small (or in the parastatistics case: when s is small), we 
nearly get Poisson statistics. In quantitative terms, the 
operator norm ||5||oo (the largest eigenvalue) measures 
well the maximal effect of statistics, in some sense the 
maximal phase space density: 



io g c s (/)-iogc (/)| < llsiliftMllsiL) 

with /^-l- 10 ^^- 



(42) 



F. For comparison: coherent states 

Quasi-free states (which are associated to e.g. ther- 
mal light) are very different from coherent states, which 
are commonly used in describing laser beams. In order 
to make the distinction clear, we derive here the char- 
acteristic function and factorial moments for the latter 
case. The (Bose) coherent states are given by the non- 
normalized vectors 



e 

N=0 



1 



(43) 



in the Bose Fock space T+(H). The characteristic func- 
tion ( 15 ) is now 



C(f) = (e*\T + (F[e*f})e*)/\\e^ 
= cxp(<j>\F{e lf - 1)4>). 



(44) 



This is precisely the characteristic function Q of a Pois- 
son random field with intensity measure 



fi(dx) = {<j>\F{dx)(j>). 



(45) 



It is remarkable that this holds for any second quantized 
observable. The factorial moments are now simply prod- 
ucts of p,: 



m k (dx x ■ ■ ■ dx k ) = fJ,(dxi) ■ ■ ■ p(dx k ). 



(46) 



Comparing with ( 33 1 , we see that a quasi-free state with 
the same first moment, i.e., with a — \cj))(<f)\, has larger 
second moment. In particular, thermal light has more 
variance than coherent light of same intensity. 

From ( 46 1 we also immediately see that for any coher- 



ent state and any counting observable, all the correlation 
functions (13) are constant, jWfii,..,,^) = 1, as ex- 
pected m jnf. 



IV. STATIONARY LIMIT OF A PARTICLE 
SOURCE 

The kind of limit we have to take here is clear from the 
case of the Boltzmann statistical model (25). We as- 
sumed there that pis a finite measure, so the total num- 
ber of particles had finite expectation. But we can also 
take p{di) — podt as a multiple of Lebesgue measure, 
where po is the emission rate, i.e., 



a = Ho / dt cr(t), 



(47) 



provided the integral is a bounded operator and, for fi- 
nite intervals S, tr ^Jo~F{S) \fo has finite trace. One can 
take this as a motivation for looking integrated trace class 
operators (47) also in the case of Bose/Fermi statistics. 
However, it is physically more realistic to have a descrip- 
tion of beam generation which is consistent with statistics 
from the outset. 



For building a stationary source model it is best to 
include the particle generation in the dynamics. In this 
way one can consider sources operating continuously for 
an arbitrarily long time. The intuitive idea is that after 
being activated at time t — 0, the source creates particles 
with fixed initial wave function </>, one after another, each 
particle subsequently evolving according to some single 
particle Hamiltonian H with direct integral decomposi- 
tion as discussed in Sect . II C| Formally, this is expressed 
(see, e.g., the review [22], and [23]) by evolving the many 
particle state p t according to the master equation 



dt 



Pi 



-i[dT(H),p t ] + \{2a^p t a^ - a^a^pt - p t a^), 

(48) 

with the initial condition that pg is the vacuum state. 
Here and in the rest of this section we consider only the 
Bosonic case. Then A > quantifies the strength of 
the source, and dT(H) is the many particle Hamiltonian 
corresponding to H. The time evolution t i— > p t is a 
quasi-free semigroup [55], so the Fock space state p t is 
quasi-free for each t > 0. This reduces the dynamics to 
a one-particle problem, which has the solution 



a (*) = 2A / ds e sT *\(t>){ct>\ 
Jo 



(49) 



where = —iH + X\(j>)(cj>\ is the generator of a strongly 
continuous semigroup. 



The integral 0^(00) is the analogue of (47). To com- 
pute it, we can formally solve the function /3^(t) = 
e(t)(e* T *<^}, in terms of 7v ,(t) = Q{t){e- UH , 
where is the Heaviside step function. In fact, we get 
$ip = S^^, where h denotes the inverse of the Fourier 
transform given in p4|, and 



(50) 



S*(E) = (1 - \%(E)Y 



This gives (V>|3>(co)V) = Xn- 1 J^JS^E^E^ 2 dE, 
corresponding to a quasi-free state which is invariant with 
respect to the evolution ( 48 ) but not with respect to the 
free Hamiltonian H . 

In order to get a state which is stationary for 
the free evolution, we now take another limit a = 
lim s _ ! . +00 e lHs a<j,(cc)e~ lHs . In the spirit of scattering 
theory, this amounts to translating the state back in 
time with the free evolution, after having let it evolve 
a long time according to the particle generating semi- 
group evolution. Since j e -iHs^(t) = 7^,(4 — s), we have 
e lsE %^Hs^(E) -> 2n((j)(E)\iP(E)) as s -> +00, where 
4>(E) is the wave function <j) in the H (energy) represen- 
tation (20). Here we have used the fact that 



itH 



HE 



(4(EM(E)) ni 



Thus the final stationary limit is given by 



d = An\ dE\S^{E)\ £ \(j )E ){4> E \ 



(51) 



(52) 
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In the ID case with the free Hamiltonian, one can al- 
ternatively take a limit of moving the source to — oo in 
space, which results in a similar expression but with the 
projection onto positive momenta. Even in this case the 
denominator, which reflects the phase space density at 
the source, still depends on the negative momentum com- 
ponents of (f>. 

There are, of course, some assumptions needed to make 
the above derivation valid. Physically, we expect [23J that 
the free evolution H should be fast enough compared 
to the strength of the source so that the particles do 
not accumulate, or even condense, near the source, but 
move away as new ones are created. Mathematically, the 
relevant assumptions can be expressed as follows: (i) The 
wave function (f> is bounded in the energy representation, 
and (ii) the operator of multiplication by S^(E) defines 
a bounded operator which keeps the Hardy class H 2 ~ 
invariant. The assumption (i ) en sures that ||7v>||2/||V'll 
is uniformly bounded due to (51); in particular, the L 2 
Fourier transform 7^, belongs to the Hardy class H 2 ~ , 
and so can be extended to an analytic function in the 
open lower half plane. Then (ii) implies that S^-y^p is 
square integrable, and has support on [0,oo); hence S^, 
is well defined, the formal relation = makes 



sense, and (52) is bounded, the limits existing in the 



weak operator topology. 

The assumption (ii) can be replaced with stronger but 
more easily verifiable versions: for instance, if 7^ is (ab- 
solutely) integrable, and y^(E) 7^ A -1 for all E in the 
closed lower half plane, then (ii) holds. Even stronger 
condition [33] is A J °° |7^(i)| dt < 1. 

The assumption (ii) implies, in particular, that 
/ I S<j,(E) 1 2 1| (j) E \\ 2 dE < 00. In the next section we will 
see in a more general context that this condition en- 
sures tr y/aF[f]y/a < 00 for any arrival time observable 
F as in Sect. |II C] and any / compactly supported in the 
time direction. By Lemma[T] the local counting statistics 
of T + F is therefore well-defined for the stationary state 
(52), and can be obtained from (36). 



In order to complete the discussion on this sta- 
tionary limit, we have to check that the counting 
statistics for the limit state ( 52 ) can be approximated 



by measuring the actual finite beam emitted by the 
source. It is clear by construction that we can write 
a s , t = e lHs d^{t)e- lHs = W* 4 W Sjt , where (W S ^)(E) = 

^/YpKe lsE $le~ iEt ' P^-isH^dt 1 . Here W s , t converge 
strongly when we take first t — > 00 and t hen s — > 00. 
Moreover, a Syt < ||S , 0|| 2 (To, where cto is (152]) without 
\S l f > (E)\ 2 . Hence, it follows from Lemma [2] that 



lim lim C s , t (f) =dct(l 

s— yoo t—^oo 



^F[e i} - l]V(?)-\ (53) 



where C Sjt is the characteristic function (36) correspond- 
ing to the operator er Sjt . 



BEAMS AND RATES 



Beams 



The following general scheme emerges from the above. 
We consider systems with Hamiltonian H, and decom- 
pose the Hilbert space into a direct integral H = 
J dEHs over the spectrum of H. Stationary beams 
are described by an extension of quasi-free states, given 
in terms of a one-particle operator 



0= I dE d(E), 



(54) 



commuting with the Hamiltonian. Here <?(£') is a posi- 
tive trace class operator in the multiplicity space at 
E, which depends on the details of the source. The basic 
normalization condition for these operators is that 
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i [dE tr a(E) < 00. 
2ir J 



(55) 



We emphasize that no such operator is trace class. In- 
deed, a general trace class operator T would have an 
integral kernel T(E, £'):%^W E so that (Tip)(E) = 
JdE' T(E,E')ij)(E ). The trace of such an operator is 
trT = JdE trT(E.E). In contrast, the operator Q 
has the formal integral kernel T(E, E') = a(E)5(E-E T ), 
which is singular on the diagonal. Such operators do 
arise from integration of trace class operators over time 
in the sense of (47): The time evolved operator U^TUt 
has integral kernel T(E,E')exp(it(E — E')). Integrat- 
ing this with respect to time we get the kernel 2ir5(E — 
E')T(E, E). Hence the trace class condition for T turns 
into (55| for the integral JdtU^TU t . 



Thedirect integral form ( [54] ) or, in other words, the 
elimination of off-energy-diagonal terms in the kernel 
for a is, of course, just the consequence of stationarity 
[a, H] = 0, and leads to a major simplification in the 
computation of expectation values and rates. 



B. Rates 

We now want to combine the stationary sources given 
by a of the form ( 54 ) , with a general counting observable 
T S F, i.e., the second quantization of a general arrival 



time observable F, as discussed in Sect. II C The full local 
counting statistics is then contained in the characteristic 
function C(f) restricted to test functions / which have 
compact support in the time direction. Technically, this 
approach is based on the discussion in Sect. HI D[ which 
guarantees the existence of C(f) in (36) for any / with 
f(x,t) — for t £ [ti,t2], once we have 



trVaF[x]Vd < 00, 



(56) 



where x(x,t) — 1 for t G [ti,t 2 ] and x( x ^) = other- 
wise. One of the consequences of this approach is that 
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we can always get a justification of the formula starting 
from finitely extended beams (trace class a) and going 
to a stationary limit in (36). There are many ways to 



But 



do such a limit, which corresponds to the many ways a 
beam which looks basically stationary during a fixed time 
interval could begin in the distant past and end in the 
far future. The formula (36) thus captures the essence of 



what we mean by "stationary beams" . 

The rest of this subsection will be devoted to substan- 



tiating the above claim (56) for a of the form (54), any 
time-covariant F, and any t\ < t 2 - We do this by showing 
the more general estimate 



/dF[fY 



< ||/||oo7l*a-ti|, 



(57) 



for any bounded complex (measurable) test function / 
such that f(t,x) = for t ^ [ti^L where 7 is the rate 
constant from ( 55 ). First we apply the triangle inequality 



to a sum of positive operators, and use that on positive 
elements the trace norm is just the trace. That is for 
positive operators F a and f a £ C we have 



< max f a I 

1 m 



E trF < 



(58) 

Hence, for step functions / = Y^ n f a Xai with X — 
SqXqi the left hand side o f J57| is bounded by 
H/lloo tr VaF[x]Vd- To prove (|57|), it is therefore suf- 
ficient to show that 



trV?F[x]V?< 7|i a -ii|. 



(59) 



We do this by expressing F by its dilation F = V*FV. 
Since F[x] is a projection, the trace we need to compute 
is tr W*W with 



W = F[ X ]V\ 



(60) 



Since tr W* W — tr WW* for any operator (where both 
sides might still be infinite) , we now compute tr WW* — 
F[x]VaV*F[x]. Note that 



VdV* = j dEV E d(E)V E 



(61) 



commutes with the energy. Therefore, in the time do- 
main it acts as a convolution operator: 

(VdV*$)(t) = [ ds S(t - s)V(s) 



with 



5(i) 



1 

2^ 



dE e ltE V E a(E)V E . (62) 



Here the integral defining 5 is convergent in trace 
norm by assumption (55 1, and by the Ricmann-Lcbesgue 



Lemma t 1— > S(t) £ &QC) is a continuous function van- 
ishing at infinity. Due to the continuity we can evaluate 
the trace as an integral on the diagonal of the kernel 
K{t, s) = S(t - s) Thm. 3.9.]: 



tvF[ X }VdV*F[ X ] = / dttv K S{0) 



(63) 



tr* 5(0) 



^ (dE trd{E)V E V E < 1 . 
2ir J 



(64) 



This completes the proof of the estimate ( 57 ) . 

According to Sect. |IIID| in order to approximate the 
counting statistics of a stationary state (described by a 
of the form ( 54 )) by the counting statistics of finite-beam 
ones (with trace class cr a ), it is sufficient to find one ma- 
jorizing (Tq of the form (54), with a a < ctq, and have 



y/Oa — > v a at least weakly. This will give uniform con- 
vergence of characteristic functions for the local counting 
statistics. 



C. Moments 

The moments are best expressed in terms of the oper- 
ator valued function S(t), defined in (62). This depends 
both on the source via a and on the observable chosen, 
via V. It also contains the required Fourier transforma- 
tions, so all moments are immediately expressed in the 
time domain. The idea is to reduce the factorial mo- 



ments to the measures from ( 34 ) , rewritten by using the 
dilation: 

Hl{dxi ■ ■ ■ dx e ) = tr (VaV*F(dxi) ■ ■ ■ VdV*F(dx e ) 
The operator under the trace has integral kernel 
K(t ,t e ) = Jdtx- ■■dt i - l S{t Q - ti)F{dx-L) x 

xS(ti -t 2 )---S(t e -i -tt)F(dx t ) 

Then the required trace is J dt tr K(t, t), where the trace 
in the integrand is over K,. In all these expressions the 
argument of the measure F is still the combination of 
time and arrival location. But noting that with respect to 
time F is just a multiplication operator, so for f(t,x) — 
h(t)g(x) we have (F[hg]$)(t) = h(t)G[g]$(t). Here and 
in the following x stands only for the arrival location. If 
we now take h and g as the indicator functions of small 
sets dt C K and dx C X, we find that 

[ie(dti dxi, ■ ■ ■ ,dti dxi) = dt\ ■ ■ ■ dtt x (65) 

x ix(s(t e - t 1 )G{dx 1 )S{h - t 2 ) ■ • • G(dxifj 

For the first (factorial) moment we thus get 



rhi(dtdx) = dt tr S{0)G{dx) 



(66) 



Hence 5(0) serves as a "density matrix" for count rates. 
It is normalized to the total particle rate 7. For a dis- 
crete family of counters with POVM elements G x > 0, 
^2 X G x = 1 the arrival rate at counter x is j x — 
tr5(0)G x . 

The second moment has a density depending only on 
the time difference r = t 2 — t\. For discrete counters we 
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FIG. 1. (Color online) Second order normalized correla- 
tion function g (t). Upper figure: Lorentzian, lower fig- 
ure: Gaussian initial spectral density. Boltzmann model 



(thick solid line), quasi- free state given by (52 I with rates 
7 = 0.5, 1.0, 1.5 (bottom to top curves). 



also give the form of the normalized correlation function 
9 



(2). 



fh\{dt\ dxi 7 dt 2 dx 2 ) 
M T {dx\, dx2) 



dhdt2M tl -t 2 (dxi,dx2) (67) 
(trS , (0)G(dxi))(trS'(0)G(da; 2 )) 
+str S(T)*G(dx 1 )S(T)G(dx 2 ) 
9%}(r) = 1 + — tTS(T)*G x S(T)G y 

Ixiy 

Here we used the symmetry S(—t) = S(t)* to write the 
expression in a more obviously positive form. 

It is not very enlightening to write down the higher mo- 
ments. The third moment (35 ) will contain contributions 
trS(t 3 - h)G x S{t x - t 2 )GyS[t 2 - t 3 )G z . 



VI. EXAMPLES 

A. Second order correlation function 

In the case of the free Hamiltonian in one dimension, 
we have He — C 2 . Taking Vb = 1 for E > 0, the source 



( 52 1 , and one detector G corresponding to the detection 



of right going (positive momenta) particles, we get 
7 = Or)- 1 ^), g^(r) = 1 + | X (t)| 2 x(0)- 2 , 

where x (r) = 4ttA / °° dEe- irB \<f> E ,+ \ 2 I 1 ~ X K E )\~ 2 , 
and <j)E + is the positive momentum component of 4>e- 
In Fig. jl] g^(r) is shown for a Lorentzian I^b^+I 2 ~ 
{a/2)/((E-E a ) 2 + (a/2) 2 )/Tr and a Gaussian |0 B ,+ | 2 ~ 
exp (— (E — E ) 2 / (2a 2 )) /(y/2lra), with no negative mo- 
mentum components. In the case of the Lorentzian the 
correlation function can be approximated for En 3> a by 
_l_ e -r(a-2\) ^ gee solid lines in Fig. Ill), i.e., as 



an exponential modified by the intensity parameter A, as 
Bose statistical effects become more relevant. 



B. Plane wave beams 

The energy density of the beam is tra(E), which 
clearly needs to be integrable. In the Fermi case the 
constraint a < 1 excludes singularities in this density. 
However, in the Bose case, we can also consider singu- 
lar distributions. Let us assume a one-dimensional, free 
Hamiltonian and we are restrict to only positive momenta 
to simplify the notation. In this case the beam state is 
given by 



dEa(E)\E)(E\ . 



where \E\ are the generalized energy eigenvectors. Con- 
sider a sequence of functions a n with a n (E) K ^(£;_ 
Eq) =: a(E). The characteristic function in this case is 

C(f) = (1 - K(E \F[e*f - l}^))- 1 . 

The rate is then given by 7q = k(E \F [f]\E ). We can 
get the number distributions p n for a measurement result 
in an interval Y from the characteristic function, see ( 10 1. 
This characteristic function for for a(E) = k8(E — Eq) is 

CiMv) = ^-j^ 1 



l-{e"-l)K{Eo\F{Y)\E ) 



OO 

E 

n=0 



(i + g) r 



(68) 



where q = k(E \F(I)\E ) . By comparing (|68| and (10 1, 



we get the number distribution for a detection in the in- 
terval Y, namely pQ. n — ■ Fig. j2j shows examples 
of the number distribution for an arrival-time measure- 
ment with q = \/T0. 

It is illustrative to compare it with the number statis- 
tics of a coherent beam given by (44) with <f>(E) = 
•Jk5(E — Eq). The characteristic function is then C(f) = 
exp (K(Eo\\F[e i f — l]|i?o)) and from this we get the rate 
7c = K(E \F[f]\E ) = 7q and a Poisson number dis- 
tribution pc.n — ^jl™ 6 ~ q which is also shown in Fig. [2] 

Let us look at a Kijowski's arrival time measurement 
in more detail. In that case, we get the rate jk = 5^ 
and qx — IkKY) where l(Y) = J Y dt. For a quasi-free 
beam we get for the probability for no detection in an 
interval [ti,t 2 ] is Pq,o(*i,*2) = 1/(1 + 7if(*2 - and 
for a coherent beam pc,o = e~ 1K ^ t2,t '^ lq . The waiting 



time calculated by ( 12 ) is now for a quasi-free beam 



wq(t) 



2jk 



(1 + Jkt) 3 



(69) 



and for a coherent beam 



9 



(2) 



w c {T)= lK e^ KT (70) 
which is -as expected- an exponential distribution. 
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FIG. 2. (Color online) Particle number distribution p n with 
q = VTO, quasi-free beam PQ, n (blue boxes), coherent beam 
Pc,n (red circles). 
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